RIKEN-TH-58 



Low energy states of a semifiexible polymer chain with attraction and the whip-toroid 

transitions 

Y. Ishimoto 

Theoretical Physics Laboratory, RIKEN, Wako 351-0198, Japan 
N. Kikuchi 

Institut fur Physik, Johannes Gutenberg- Universitat Mainz, Staudinger Weg 7, D-55099 Mainz, Germany 

(Dated: February 2, 2008) 

We establish a general model for the whip-toroid transitions of a semifiexible homopolymer chain 
using the path integral method and the 0(3) nonlinear sigma model on a line segment with the local 
inextensibility constraint. We exactly solve the energy levels of classical solutions, and show that 
some of its classical configurations exhibit toroidal forms, and the system has phase transitions from 
a whip to toroidal states with a conformation parameter c = ^ (^) • We also discuss the stability 
of the toroid states and propose the low-energy effective Green function. Finally, with the finite size 
effect on the toroid states, predicted toroidal properties are successfully compared to experimental 
results of DNA condensation. 
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I. INTRODUCTION 

In nature, biological macromolecules are often found 
in collapsed states 0, 0, Q • Proteins take unique three- 
dimensional conformations in the lowest energy state (na- 
tive state), which is of great importance in its function- 
ality 0. DNA in living cells is often packaged tightly, 
for instance, inside phage capsids. Recent advances in 
experimental techniques mean it is now possible to study 
the conformational properties of biopolymers at single 
molecular level 

El El Hi 01- 

As well as its biochemical, 
medical, and industrial importance, (bio-)pol yme rs have 
drawn much attention 0, 0, IM S 13, IM. S [lO, HH E3, E3, 

MM 

To increase our understanding of their physical prop- 
erties, a flexible homopolymer chain in a dilute solu- 
tion, as the simp lest model, has been heavily investigated 
QIaSElHEHIlllill. When the temperature 
is lowered, or the solvent quality is changed from good 
to poor, the resulting effective attractive interactions be- 
tween monomers can cause the polymer to undergo a coil- 
globule transition (collapse transition) from an extended 
coil to a compact globule state [H E. 01. Both equilibrium 
0, II II EE 03 and dynamical [Sll! HI EI properties 
of the coil-globule transition of the flexible chain are now 
well understood. 

However, many biological macromolecules such as 
DNA, F-actin, and collagen show large persistence 
lengths and are classified as semifiexible chains p H23 . l23| . 
For instance, double stranded DNA in aqueous solution, 
mostly with segment diameter a ~ 2 nm, has the persis- 
tence length / ~ 50 ~ 60 nm. Therefore natural DNAs be- 
have as semifiexible chains when their contour lengths are 
several orders longer than I 0, 122L |23| . In such cases in a 
poor solvent condition, the balance between the bending 
stiffness and surface free energies induces toroidal confor- 
mation rather tha n sp herical glob ule of a flexible chain 
13 El El IHEi El IS HSi. In fact, when we 



put condensing agents as multivalent cations into DNA 
solution, it can cause DNA to undergo the condensation 
from a worm-like chain (whip or coil) to toroidal states 

cannula. 

Towards the understanding of the "whip (or coil)-toroid 
transition" of a semifiexible homopolymer chain, or of a 
DNA chain, many experimental and theoretical works 
have been done, in particular, in a poor solvent condi- 

tioqJjj-PU^Jjl, 1 !!; jlBlaaiallalEllalEl 

llallallaMliinilMHlil. Extensive results from 
experiments showed that collapsed DNA exists in toroid, 
rod, sphere and spool-like pha ses with the toroid being 
the most probable [M El El El EI Simulations us- 
ing Monte Carlo, Langevin approaches or Gaussian varia- 
tional method, calculated phase diagram for the semifiex- 
ible chain in a poor solvent HmiSIHEllMElEH- 
In theoretical works, existing phenomenological models 
balance the bendin g an d surface free ener gies to estimate 
toroidal properties jH |H |H m HI Elll IH H H| . 
It becomes increasingly probable that toroid is the stable 
lowest energy state — the ground state. 

We note, however, that the theoretical aspects of the 
works assume a priori toroidal geometry as the stable 
lowest energy state with no theoretical proof |2l . More- 
over, compared to the theory of coil-globule transition 
of a flexible chain 0, H H 0, 0, which are well de- 
scribed by Gaussian approximation and field theoretical 
formalism [ill IT^. there is no simple "microscopic" the- 
ory, which contains the salient physics to demonstrate 
the whip-toroid transition of the semifiexible polymer. 

Difficulties in formulating theory results specifically 
from the "local inextensibility constraint" of the semi- 
flexible chain, which makes the theory non-Gaussian , 
and also from the "non-local nature" of the attractive 
interaction along the polymer chain, which makes the 
theory analytically intractable. As a result, even for the 
simplest semifiexible chain model without attraction, i.e. 
the Hamiltonian J3J|, only a few equilibrium properties 
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are analytically tractable such as the mean square end- 
to-end distance (R 2 ) of a free chain [53, and that of 
a semiflexible chain confined to a spherical surface ^3] . 

To overcome these problems, we propose a microscopic 
model to describe the whip-toroid transitions of a semi- 
flexible homopolymer chain at low energy — at low tem- 
perature or at large persistence length. To explore the 
equilibrium distribution (Green function) of a semiflexi- 
ble chain, the path integral formulation is applied rather 
conventionally. Note that a semiflexible homopolymer 
chain in equilibrium at low energy satisfies the local in- 
extensibility constraint. Also, if the chain satisfies the 
local inextensibility constraint, its Hamiltonian becomes 
equivalent to the 0(3) nonlinear sigma model on a line 
segment. Therefore, a semiflexible homopolymer chain 
at low energy can be formulated in the path integral of 
the 0(3) nonlinear sigma model on a line segment. It is 
the first time that the local inextensibility constraint and 
the non-local attraction in the path integral are employed 
together and are solved clearly. Exploring it in detail, we 
find the toroid states as the ground state and the whip- 
toroid transitions of the semiflexible chain at low energy, 
which can also be found in our preprint |48j . We then 
discuss and test the stability of the toroidal solutions, 
and propose the low-energy effective Green function. We 
show, in final sections, that our predictions on toroidal 
properties are in sufficiently quantitative agreement with 
the experiments [T^.[l^|. 

The paper is organized as follows. In sections II and 
III, a semiflexible polymer chain with a delta-function 
attractive potential is formulated in the path integral 
method. We then deduce 0(3) nonlinear sigma model on 
a line segment with the local inextensibility constraint. 
In section IV, we derive the classical equations of mo- 
tion for the nonlinear sigma model action, and solve 
them explicitly. We also prove that our solutions rep- 
resent the general solutions of the equations. The pre- 
cise microscopic Hamiltonian, or the energy levels, are 
obtained from the solutions, and the conditions for the 
stable toroids are given. We also investigate the phase 
transitions in the presence of the attractive interactions. 
Section V is devoted to the stability of the toroidal states 
under the 'quantum' fluctuations away from classical so- 
lutions. We also construct the low-energy effective Green 
function from those of the whip and toroid states using 
perturbation theory. In section VI, the finite size effect is 
introduced and the theory is mapped onto physical sys- 
tems. Assuming the hexagonally packed cross sections 
and van der Waals interactions, we show that our micro- 
scopic model does fit well quantitatively with a macro- 
scopic property of the toroids — the mean toroidal radius 
in the experiments 0, 0] . In the final section, our con- 
clusion summarizes the paper and discussions are given 
with respect to the literature and the future prospects. 
Note the precise definition of the delta function potential 
is given in Appendix A, and the SO (3) transformations 
are described in Appendix B. 



II. POLYMER CHAIN AS A LINE SEGMENT 

In the continuum limit, the Green function (end-to- 
end distribution) of a semiflexible polymer chain with 
attractive interactions can be given by the path integral: 

G(Q,R-u u u f -L,W)=N- 1 [ r{L) 5[ffs)]e- w [ f,s '^ (1) 

Jr(a)=8,u(0)=Ui 

with the local inextensibility constraint \u\ 2 = 1 p2l.l23|. 
s is the proper time along the semiflexible polymer chain 
of total contour length L. r(s) denotes the pointing vec- 
tor at the 'time' s in our three dimensional space while 
u(s) = corresponds to the unit bond (or tangent) 



vector at s. J\f is the normalisation constant (E 
Following Freed et al. and Kleinert [rH [22 
mensionless Hamiltonian can be written by 



H[f, u, W] 



ds [H(s) + V AT {s)} 



the di- 



(2) 



where H(s) and Vat(s) are the local free Hamiltonian 
and the attractive interaction term, respectively: 



His) = - 
V AT (s) = 



d 

OS 



-W I ds'5(f(s)-f(s')). 
'o 



(3) 
(4) 



I is the persistence length and W is a positive coupling 
constant of the attractive interaction between polymer 
segments. Thermodynamic j3 — l/(fcsT) is implicitly 
included in I and W, which can be revived when we con- 
sider the thermodynamic behaviours of the system. I is 
assumed to be large enough to realise its stiffness: l^> lb, 
where lb is the bond length. Note that there has been no 
consensus about the form of attractions, but people in 
the literature agree that effective attractions derive the 
toroidal geometry [Ij|. For example, in DNA conden- 
sations, interplay between charges, salt and other un- 
settled (unknown) elements derives extraordinary short- 
range dominant effective attraction in a poor solvent con- 
dition. Therefore, we introduce the above delta-function 
potential Vat{s) for the modelling of the DNA conden- 
sation in a poor solvent condition, again as in Freed et 
al. [Hill. As you can read off from the above, Vat(s) 
takes the non-local form, since the form at s contains in- 
formation at the other points s' £ (0, s). In Vat(s), we 
omit the symbol for the absolute value \f(s) — r(s')\. (see 
Appendix A for the precise definition of the potential.) 

In what follows, we express r by the unit bond vector u 
and therefore the Hamiltonian TL{u) in terms of u. Hence, 

the Green function G ^0, R; iti, u/; L, W^j becomes a path 
integral over u with the positive coupling constant W, 
regardless of r, 



G 



V\u 



( s )\ 6 (lo ds u( s ) 



R 



-H[U,W] 



(5) 
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where we used r(L) = J ds u(s) and the Jacobian is ab- 
sorbed by N which is neglected here. The delta function 
selects out the end-to-end vector. Basic properties of the 
Green function is given below. 

Due to the local inextensibility constraint |u(s)| 2 = 
|<9r(s)| 2 = 1, the total length of the polymer chain is 
strictly L for G(R, L, W). Thus, the Green function as a 
distribution function exhibits a hard shell at \R\ = L: 



G (o, R;Ui,u f ;L,W 



for \R\ > L. 



That is 



\R\<L 



d A RG [0,R; Ui,u f ; L,W) = 1 



(6) 



(7) 



It further means that the normalisation constant is given 
by 

N = I d 3 R ( U! V[u{s)\ S frf dsu(s) - R) e- u ^' w \(8) 

J\R\<L Ju, ^ 1 



III. 0(3) NONLINEAR SIGMA MODEL ON A 
LINE SEGMENT 

When W = 0, our free dimensionless Hamiltonian is 
given solely by u field: 



H(u) 



H{r, u, W 

1 r L 

2 Jo 



0) 



ds \du(s)\' 



(9) 



with the constraint |u(s)| 2 = 1. This can be interpreted 
as the low energy limit of a linear sigma model on a 
line segment, or quantum equivalently a nonlinear sigma 
model on a line segment, rather than some constrained 
Hamiltonian system. 

In this section we consider 0(3) nonlinear sigma model 
on a line segment for the path integral formulation of 
the semiflexible polymer chain. This is nothing but a 
quantum mechanics of a limited time s € [0, L] with a 
constraint. The constraint \u\ 2 = 1 restricts the value 
of u on a unit sphere S 2 . This can be transformed into 
u 2 — 1 — u\ - u\. 

Substituting this into eq.@ gives 



S[ui,u 2 ] 



ds [G ij ] dui(s)duj(s) (10) 



where the metric G lJ on the unit sphere in three dimen- 
sional w-space 



G K 



[u\,u 2 \ 




(11) 



This is called the nonlinear sigma model since the action 
is 0(3) symmetric but some of its transformations are 



realised nonlinearly on this {ui} basis. It is also equiva- 
lent to the classical Heisenberg model with a constraint 

of unit length spins Si = 1 in the continuum limit |50j . 

The action can also be expressed in the polar coordi- 
nate: 



ui = r u sin 6 U cos tp u 
u 2 = r u sin U sin ip u <4> 
"3 = r u cos9 u 



r u = \u\ 

6 U = arccos^ ,(12) 
= arctan — 



'•P 



S[6 w <Pu] = o ds [(d9 u ) 2 +sin 2 e u (d^ u ) 2 



ds [G u ] d6i(s)d0i{s) 



(13) 



where (9 1: 9 2 ) = (9 u ,ip u ), and the metric G y is given by 
the diagonal matrix: 



G* j [0i,0 2 ] 



1 
sin 2 6»i 



(14) 



This is essentially the same as G u [ui\ since both are met- 
rics on the same sphere S 2 . The 50(3) transformations 
of the polar coordinates (9 u ,ip u ) can be expressed by 
three infinitesimal parameters gi (see Appendix B): 



59 u = gi sin cf> u - g 2 cos <j> u , 
8(j>u = cot 9 U (gi cos 4> u + g 2 sin ( 



93- 



(15) 



The canonical quantisation of the action 113|) is suit- 
able for the investigation of the local nature of the sys- 
tem, but not for its global nature such as toroidal con- 
formations. Therefore, we focus on the classical solu- 
tions of the action i|13fl and consider the quantum fluctu- 
ations around the classical solutions using the path inte- 
gral method. Integrating the action (|13f) by parts gives 



I f L 

S[9 U , ip u ] =-- J ds [9 u d 2 9 u + (p u (d o sin 2 0„ o d) Lp y 



where 



Surface 



stands for the 



(16) 



pings and the surface term 

1 L 



composition of the map- 

L 



Surface 



1ud9 u 



(i 



The surface term might be neglected 



by taking the north pole of the polar coordinates 
(9 U (0), (p u (0)) = (0,0) and considering the static solu- 
tions. By setting the north pole, half of the surface term 
vanishes. Given that we have the static solutions, i.e. 
u(s) ~ (u) the surface term contribution becomes much 
smaller compared to the bulk term ~0 Of) where lb is 
the constant bond length. Minimizing the action i|13|) 
in terms of 9 U and ip u yields the classical equations of 
motion: 



29,, 



= 



"a 2 + 2(<90 n )cot0„<9] ip u = 0. 



(17) 
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IV. CLASSICAL SOLUTIONS AND THE 
WHIP-TOROID TRANSITION 

Our aim in this section is to explore classical solutions 
of eq. lfP7|l and to study the lowest energy states and the 
whip-toroid phase transition in the presence of attractive 
interactions. 



A. Classical solutions 

Consider classical solutions of ca. H17|) with a trial so- 
lution 9 U — 0. The first equation of (|17|l leads to 
sin 29 u ((fiu) 2 = . Thus, the solution is either 9 U = , \ , 7r 
or ip u = 0. The solutions 9 U = 0, n or ip u = with 9 U = 
are equivalent to having a constant u. Accordingly, clas- 
sical solutions reduce to 9 U = ? or ti = const. When 
we substitute 9 U = £ into the second equation of motion 
(fP7|l . we obtain d 2 <p u — 0. Therefore, we have the two 
classical solutions 

u(s) = const, 
or 

IT 

9 = — and ip u — as + 6, (18) 

where a, b are constants. Note that the second classical 
solution of eq. I|18|l is the uniform motion of a free particle 
on the sphere (see Fig0. 

By symmetry argument, we state that the solutions 
(|18|l represent all the classical solutions. That is either a 
constant u(s) (rod solution) or a rotation at a constant 
speed along a great circle on the S 2 (toroid solution). 




FIG. 1: Classical solutions of ea. l|17|l : (a) constant u, (b) a 
path along a great circle on S 2 . 

Proof) 

The theory has 0(3) D 50(3) global symmetry. Ac- 
cordingly, one can take any initial value of u(0) for a 
classical solution. In other words, one may set u(0) to 



be the north pole for the representatives of the classi- 
cal solutions, using two degrees of freedom of SO(3) ro- 
tations. In addition, by 50(2) local rotation symme- 
try or by one residual degree of freedom of 50(3), we 
can freely set the orientation of du(0). For example, 

(^9 u (0),ip u (0)j = (0,a) or (a,0). So, one may set the 

initial values as 

u(0) = (0„(O)>¥>»(O))=(~,O), 
du(0) = (fl„(0),6,(0))=(0,o), (19) 

with a condition a > 0. The non-negative real constant a 
turns out to be the only degree of freedom that represents 
all the classical solutions. 

Substituting these initial values to the equations of mo- 
tion H17|). we obtain at s = 

d 2 9 u = 0, 0V = 0. (20) 

So, an infinitesimal change e of the variable s yields 

(9 u (e), p„(e)) - (|, eo) , (V u (e), £,(e)) - (0, a). (21) 

As one can see in the equations of motion and in the 
above, so long as 9 u (s) = J, ea. ij2U|) holds at any s. 
Hence, the initial conditions leads to the pair of condi- 
tions, d9 u (s) = and 9 U (0) — -|. In other words, the 
pair of the conditions exhaust the representatives of the 
classical solutions. Thus, found solutions may well be 
regarded as the general solutions. (Q.E.D.) 

Note that the solutions l|18|) can be regarded as 'topo- 
logical' solutions in a sense that they are solitonic solu- 
tions. 



B. Non-local attractive interactions as a 
topological term 

Now we consider the attractive interaction term Q . It 
is difficult to interpret it in the context of quantum the- 
ory due to its non-local nature along the polymer chain. 
However, we can solve them with our classical solutions 
(JTSJ. Let us rewrite eq.(QJ with 

r(s)-r{s')= [ dtu(t)-( dtu(t) = f dtu(t), 

JO JO J s' 

(22) 

that is, 

Vat(s) = ~ W J o ds ' 6 dtu(t)j . (23) 

Hence the problem is now reduced to the one in the u 
space: finding non-zero values of 8 (Jj 5 , dtu(t)) with the 
classical solutions (|18fl . That is to find u(s') for a given 
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s, which satisfies \J^, dtu\ — 0. Note that, exactly speak- 
ing, the integration over s' is from to s — e with an 
infinitesimal positive constant e (see Appendix A). Thus, 
we exclude the s' — s case in the following. 

In the polar coordinates l)12[i . this is expressed by 



dt sin 9 U cos ip u = 0, 
(it sin U sin <p u = 0, 
cos^u = 0. 



(24) 



The first classical solution (u — const.) does not satisfy 
these equations and thus derives no attractive interac- 
tions. If we substitute the second classical solution of 
eq.lQSJ into eas. (|2*4*jl . we have cos#„(s) = 0, 

/ dt cos(at+b) = - (sin(as + 6) — sin(as'+6)) = 0, 
Js> a 



dt sm(at+b) = - (cos(as'+6)-cos(as+6)) =0.(25) 



Hence we have solutions: s — s' = 2nir/a > 0, n G Z. 
Without any loss of generality, we assume a > and 
n G Z + . Introducing N(s) = [as/27r] by Gauss' symbol 
EH, we obtain 



dt u{t) 



S — 1-K J a 



dtu(t) 

— Art /a 

dtu(t)=0. (26) 



Therefore, the attractive potential is given by 
Vat(s) = -W-N(s). 



(27) 



Note that N(L) represents the winding number of the 
classical solution (fl%|) along a great circle of S 2 (see 
FigD- Finally, an integration over s yields the dimen- 
sionless Hamiltonian with our classical solutions: 



H [u,W} = / dsH(s)+ / dsV AT {s) 
Jo Jo 



LI , 
= — a 2 -W 
2 

LI 



O N W- 1 o / T 



fc=l 



y a 2 ~iyi.7V(i){l-^(A(L) + l)} 



(28) 



The first term denotes the bending energy, and the sec- 
ond and the third terms are thought of as 'topological' 
terms from the winding number. When the chain of con- 
tour length L winds N(L) times we have the N(L) circles 
of each length ^ and the rest (L - 2zLjV(£)). The sec- 
ond and third terms in the second line of ea. l|28|l result 
from the former and the latter respectively. 



C. The toroid and whip states 

The non-zero winding number of the classical solution 
in the u space means that the polymer chain winds in 
the r space as well. That is, when a > configurations 
around the second classical solution (|18jl start forming a 
toroidal shape since 



i {sin(as + b) — sin(6)} 
f{s)— \ — i {cos(as + b) — cos(6)} 
const. 



(29) 



and stabilise itself by attracting neighbouring segments. 
We call such classical solutions the "toroid states." 
Whenever a increases and passes through the point 
for n G Z_|_, another toroid state appears with the in- 
creased winding number n. Note that the radius of the 
toroid state is given by - (see Kig- EJ) - When < a < ^5, 



Whip states 



Toroid states 





N= 1 



N = 2 



N=0 

FIG. 2: The whip (N = 0) and toroid states (N > 1). The 
value of b is given by the initial value of the bond vector Hi. 

the chain cannot wind like the toroid states. Both ends 
of the chain are not connected to each other, thus can 
move freely as well as any other parts of the chain fluc- 
tuate. As long as the total energy of the chain does not 
exceed the bending energy of ^jt at a = 2r, they can 
whip with zero winding number. 

We call such low-energy extended coil states the "whip 
states." Although the definition includes fluctuations 
around the classical solutions, unless otherwise stated, 
we primarily refer to the classical solutions of such states, 
which are rather bowstrings than whips. 

In the next subsection we explore the exact energy lev- 
els of the whip and toroid states, and discuss the phase 
transitions between these states. 



D. Favoured vacuum and toroid-whip transition 

The dimensionless Hamiltonian of the second classical 
solution (|18fl is a function of l,L,W and a: 

Hd (a,l,L,W) = ^a 2 + —N(L) (ML) + 1)-WL- ML) -(30) 
2 a 
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This matches with the first classical solution when 
N ( L ^ = for a — is defined. Accordingly, the above 
expression is valid for all classical solutions. Note that, 
since previous works assume a priori toroidal shape, no 
one clearly derived the precise microscopic Hamiltonian. 
Thus, we are now in a position to investigate exact energy 
levels of the whip and toroid states. 

Consider first a case with L, W, and I fixed. By defi- 
nition, TL{a) = TL c i{a, I, L, W) is continuous in the entire 
region of a > and is a smooth function in each segment: 



a e 



2itN 2ir{N+l) 



L 



L 



for N £ Z> . 



(31) 



However, it is not smooth at each joint of the segments: 

§^ G Z + . Introducing a new parameter c = (^r) -jjrout 
of three existing degrees of freedom, we plot in Figyjthe 
energy levels as a function of a for different values of c, 
showing qualitative agreement with Conwell et al. for the 
condensation of 3kb DNA in various salt solutions [T3 |. 
Note that, in what follows, we call the segment H31(l the 
"N-th segment" counting from 0-th, and we also call c 
the "conformation parameter" because the parameter c 
solely determines the shape of this curve. 



Energy Level 
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The first and third cases are physically less relevant since 
they mean no (meta-)stable point in the segment. So, we 
focus on the second case. 

The condition on N for the second case turns out to 
be (see Fig. EJ) 
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Region of Integer N for Minima 

































































































































































Upper boLUK 


















Lowerr bound 



FIG. 4: The solid line is the upper bound and the dashed line 
is the lower bound of N for the minima, i.e., Nu,l{c) (and 



c {N) ) 
(c (N) 

( C U,L 



The asymptotic values of Nu,l{c) are both Nu,l(c) 
~JV). 



N L {c) < N < Nu(c) for c > 4, 

1 < N < Nu{c) for < c < 4, (32) 



where 



FIG. 3: The dependence of the energy 7i.(a) on x — aL/2n 
and c. Ti.(a) is scaled by the factor of \JWl for convenience. 

Suppose N(L) = N is fixed, the Hamiltonian i|30|) 

takes a minimum at a = a c (N) = (^f N(N + 1)) 1/3 . 
Accordingly, each segment falls into one of the following 
three cases: 

(i) When a c (N) < ^j^-, H(a) is a monotonic function 
in the segment and takes its minimum at a — ^j^- . 

(ii) when 2iJV < a ^ < 2,{N+i) ^ H ^ behaves 
quadratic in a and takes its minimum at a = a c (N). 

(iii) When 27r ( j y+ 1 ) < a c (N), Ti(a) is monotonic in the 



N L {c) 
Nu(c) 




(33) 



segment and takes its minimum at a 



27r(jV+l) 



Note that, by replacing N with N(L) = [aL/27r], one can 
read the condition on a as well. 

As one can see in Figs. |3 01 , there are apparently 
more than one (meta-)stable toroid states at most val- 
ues of c. This is because the first term of bending in 
eq. (|30|l is monotonically increasing, while the other two 
terms in ea. pOl) are decreasing but not smoothly. This 
non-smoothness and the balance between two factors lead 
to multiple local minima and potential barriers between 
them. The number of minima is roughly given by the 
width of the region for N, i.e., Njj(c) — Nl(c). For ex- 
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ample, when c > 4, 



-V, ir)-X L (c } = I + C - [\jl + - c 




> 3, (34) 

where (a)k = a(a + 1)- ■ -(a + k — 1) is the Pochhammcr 
symbol. Therefore, there are at least three minima with 
positive winding numbers greater than 1. When < 
c < 4, the condition of having more than three minima 
is c > f. To summarise, when c > | there exist at 
least three minima with positive winding numbers. It 
might be helpful to mention that, if we introduce the 
finite size effect in section VI, the number of minima 
could be reduced in some cases. 

One can plot the critical value of c where the minimum 
of the iV-th segment emerges and vanishes. The lower 
bound of the N-th segment is 



N 2 



N+l 



<N, 



(35) 



> JV + 1. 



while the upper bound is 

w (N + lf 
N 

So, when c satisfies the following inequality relation: 

cf ) <c<cW ! 



(36) 



(37) 



the iV-th segment has a minimal and (meta-)stable point. 
For example, when i < c < 4, the first segment a 6 

[if' T"! (i" e " = 1) nas a minimal point at a = a c (l). 

Now we discuss the critical points of the conformation 
parameter c at which the conformational transitions be- 
tween states may occur. When Njj(c) < 1 (i.e. c < |), 
the second condition in eg. (j32|) vanishes and thus the 
whip states only survive at low energy. In this param- 
eter region, the a — rod state will be favoured as the 
ground state with vanishing energy. Including 'quantum' 
fluctuations around a = 0, we call this phase the whip 
phase. Successively, at the critical value of c = i, the 
whip phase to whip-toroid co-existence phase transition 
would occur. On the other hand, when c > ^, there al- 
ways exists at least one (meta-) stable toroid state with 
positive winding number N(L). As c grows over i, the 
local minimum in the first segment decreases from some 
positive value. Finally, when the energy of the N = 1 
stable toroid state balances with the ground state of the 



whip state (i.e. TL c i = 0), the whip-dominant to toroid- 
dominant phase transition may occur. Such a value of c 
is 27/16. Since there is a potential barrier between the 
a = rod and the N = 1 stable toroid states, the tran- 
sition is first order. When c > 27/16, the toroid states 
will dominate the action. The energy plot (Fig.[3J) clearly 
shows that the transitions between the toroid states are 
also first order, if any, as there exist potential barriers 
between two successive minima. Further discussions on 
the phase transitions will be given in the final section. 

For later convenience, we rewrite the Hamiltonian i|3U|) 
and a c (N) in terms of c and the new variable x = 



WL 

H cl (a,l,L,W) = —H(c,x) 



cH(c,x), (38) 



where 



H(c,x) = ^- + -[x]([x) + l)-2[x]. (39) 
2c x 

Therefore, [x] = N{L) and a; c ([a;]) = = 
{c-[x]([x} + l)} 1/3 . 



STABILITY, QUANTUM FLUCTUATIONS, 
AND PERTURBATIONS 



So far we have dealt with the classical solutions, which 
are derived from the first derivative of the action. Thus, 
they may correspond to the global/local minima of the 
action in the configuration space. However, the solutions 
are not necessarily stable unless we take into account the 
attraction, since the second derivative test of the action 
with W — gives the non-positive Hessian. That is to 
say, they seem to be saddle points. 



S 2 S[d u ,^ u ] 



S 2 S[9 U ,^ U } 



501 



detH 



(lip 2 ) cos26> M , 



-I (29 u (f>u cos(26>„) + Cp u sin(26> u )) , 



d^S d v d e S 
d e d v S d 2 S 



< 



(40) 



In fact, the general whip states do not need to live in a 
flat plane in R 3 whereas the classical whip state does. 
So, the transitions between the classical and the non- 
classical whip states have the flat directions, i.e., they 
can be seamless without any change of energy. Therefore, 
the stability problem is to be treated carefully with and 
without attraction. 
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A. Stability and quantum fluctuations with 
attraction 

When the attraction is turned on, the toroid states 
with the winding number of more than two may become 
extremely (meta-)stable under the quantum fluctuations 
away from the classical solutions. It is not easy to show 
that all such second derivatives of the action give positive 
values and therefore stabilise the states, since the inter- 
action term contains a special function of the quantum 
variable u. However, there is a much easier way to see 
the stability. 

Consider any small fluctuation of a segment Sl s from 
such a state. It gives rise to an increase of the energy: 

SH(a) > W ■ 5l s . (41) 

More generally, we can write down the dimensionless 

Valley of the Hamiltonian in the configuration space (c=4) 

H(x,ls) 



where N is the local number of the overlapped segments 



N 



r a(L-Z 3 ) 



]. When l s — L, the energy becomes E 




FIG. 5: A sketch of the energy level Tt(x,l a ) with the attrac- 
tion. For ease of use, the maximum value is suppressed at 4 
and la is normalised by L. The direction along the valley is 
parameterised by a variable of the classical solution. Its per- 
pendicular direction is the quantum fluctuation l a by SO(2) 
away from the classical solutions. As one can see, l s is a flat 
direction for the whip states (x < 1). 

Hamiltonian as a function of a, or x, and l s , that is, 
the energy levels away from the toroid states. l s is now 
defined as the length l s of the polymer segment shifted 
from an end of the toroid. The shifted polymer segment 
is locally rotated by 5*0(2) transformation, for example, 
by 90 degrees at each node, keeping the local bending 
energy unchanged. 

H{a,l s ) = H(x,l s ) 

-WL[x(l- l j.)] + Wl s [x(l- l f)], (42) 

where x — (Fig|SJ). The energy loss for the infinitesi- 
mal segment Sl s is found to be 

6E = WNSl s , (43) 



-|^a 2 of the bending energy. 

Therefore, under perpendicular 'quantum' fluctuations 
away from the classical solutions, the toroid states which 
look stable in Ti c i are also generally (meta-)stable. It 
should be noted here that the fluctuation along the clas- 
sical solutions may still be flat so that the transitions 
from one toroid state to another is possible. These facts 
justify our claims on the phase transitions between toroid 
states and their existence. 

As for the toroid state with the winding number one, 
its stability depends on the value of c. The transition 
along the classical solutions is almost flat so that when 
c < h it naturally goes down to a whip state (Fig[3J). 
Therefore, it is absolutely unstable. When c > ^, one 
may state that it is meta-stable since attracted parts 
locally stabilise the state. However, the non-attracted 
parts are still free to move unless it gives an increase 
of the total energy. So, the toroid state with N = 1 is 
partially stable or metastable. Its probability is given 
by summing over such quantum fluctuations of non- 
attracted parts that give the energy similar to that of 
the classical solution. 

On the other hand, even with the attraction, the whip 
state is unstable since it is not affected by the presence 
of the attraction. Therefore, it would be meaningless to 
pick up any particular shape of the whips and estimate 
its probability. Instead, one should only estimate the 
probability of all the whip states that have the similar 
energy TL{a < 2n/L), by carefully counting the number of 
such states, or equivalently by estimating entropy. Note 
that, roughly speaking, the whip state is more probable 
than a single rod state with a = 0. 

With the above reasons, it would be more appropriate 
to state that one of the toroid states of N = [x] > 2 is the 
ground state when I is much larger than the bond length 
lb and c > 4 where the whip states become negligible. 
Although we listed the above reasons, we remind that 
there is a first order transition between the rod (a = 0) 
and the toroid state (N = 1). The potential barrier 
between them is given by ^j^, thus the transition will 
be suppressed by the factor of e _1 or smaller when L < 
2ir 2 l ~ 1 fun in the case of DNA with I ~ 50nm. 



B. Perturbation by the classical solutions 



In order to complete the theory at low energy, we con- 
struct the low-energy effective Green function G e f f from 
those of the toroid and the whip states, in perturbation 
theory. To make the function more accessible, we fix the 
persistence length I in what follows. Accordingly, G e // 
becomes a function of L and W, or equivalently, of c and 
L. In addition, c = ^ (^-) >4is assumed to ensure 
the existence of the stable toroid states at the beginning. 

Let us denote the Green function of the toroid states by 



Gt , and that of the whip states by G w . As we would like 
to sum over all toroid contributions to Gt, the end-to-end 
vector R and the initial and final bond vectors (£tj,t?/) 
will be omitted in Gt as well as in G w . Therefore, Gt 
is a function only of c = c(r) where c(r) is a function of 
the chain length r of the toroidal segment. G w is also a 
function of t: the length of the non-attracted whipping 
segment in this case. Note that, however, G w {t) does not 
depend on W since the chain segment is free by definition. 

With these specifications, the effective Green function 
can be constructed by the following perturbations: 

G eff {c,L) = G w (L) 

rL-L min 

G T (c(L)) + 2 / dr G t {c{L-t))G w (t) 
Jo 

+ / dn dr 2 G w (t 1 )G t (c(L-t 1 -t 2 ))G w (t 2 ) 

J T 1 >0,T 2 >0, 

T l + T 2 < L ~ L min 

L-L min 

dr Gt(c(L~t))G t (t) 
+2 f d n dT 2 G T (c(L - ti - t 2 )) G^n) G w (t 2 ) 
+ [ dr x dr 2 G T (c(L - n - r 2 )) G w (ti) G t (t 2 ) 

J T 1 >0,T 2 >L mi „, 



l,3>0.T2>£ min , 



ri dTi 



G„(ti) x 



G T c G T (T 2 )G,„(r3) 



+2 



T l,3 >°> T 2 > L m 
S?=l Ti<L-L r , 



JJdn 



G w (ti)G t (t 2 )G w (t 3 ) 



l,2,4>°>^3> L miT7 

Ef =1 T i <l-i mir 



,i=i 



Gt c L-Etr 



G U) (ti) x 



G T - 53 t,J J G w (t 2 ) G t (t 3 ) G w (t 4 ) 

+ •••, (44) 

where L TO , n is given by the lower bound of the conforma- 
tion parameter c: c{L m i n ) = 1/2, for the existence of a 
(meta-)stable toroid state. It reads 



L „ / I 

JWV) V w 



(45) 



The first bracket in ea. H44l) gives the contributions from 
the conformations which contain only one toroid, the sec- 
ond bracket gives the ones with two toroids, and so on. 
It should be noted that so-called 'tadpole' conformation 
appears in this low energy perturbation as the third term 



in ea . I|44|l with one toroid and one whip. Schematically, 
ca. (|4*4fl can be depicted in FigEl 



Zero toroid 



o 



■Toroid 



■Whip 



One toroid 

O O o~ 

Two toroids 

OOOOOO oo<x><x> 

Three toroids 

OOO OOOOOO 

oooooooooooo 
ooo<xx>ooo 

FIG. 6: Each term in ea. l)44|l is expressed by a product of 
the toroids (circles) and the whips (waves). The third term 
is so called the 'tadpole' conformation. Below the forth term, 
all contributions are of multi-tori. 

Roughly speaking, in our ideal toroid with zero thick- 
ness, the winding number N c of the dominant toroid state 
is proportional to c(r) which is a quadratic function in 
t. Besides, the minimum of the Hamiltonian is given at 



Therefore, 



a c (N c ): Mm(W(a)) ~ ~ -^R£) 

the ratio of the probabilities of the toroids with different 
lengths t and L can be estimated by 



Gt(c(t)) 
G t (c(L)) 



■kW 2 ( r \3 
? 41 \ 2tt J 

3 41 \ 2tt J 



lit 3 - 3 ) 



(46) 



Since we assume that L is large enough as L 3 » 2 ^i l , 
the toroid states with smaller contour length are highly 
suppressed by the above factor (From c>4, we have 

L 3 > (2 5 tt 2 /IF)^ (2j^iV If W~0(1) and />!, we ob- 



tain L 3 > 

For example, this condition holds for the cases of DNA 
considered in the next section. Hence, the above pertur- 
bation can be justified for a large value of c(L). Note, 
however, that the statistical weight for each path, or each 
conformation, is not specified here, nor for the whip state 
and the normalisation factor. Moreover, we have not 
counted other conformations such as conventional 'tad- 
poles' with overlaying whips. Therefore, and unfortu- 
nately, we do not go into the precise comparison of the 
above terms. Instead, we will make some remarks on 
these issues in the final section. 



VI. COMPARISON WITH EXPERIMENTS 

In previous sections, we have found the toroid states as 
the classical solutions and found that some of them are 
(met a-) stable, one of which becomes the ground state at 
large c. As mentioned, it is pointless to ask Ui,Uf, and 
R(L). One of the most physically meaningful observable 
is the radius of the toroid. 
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For large c, the ground state — the dominant toroid 
state of the winding number jV c can be estimated by the 
inequality relation (|37[1 of c: 



< c < c 



(ATc 



By its inverse relation, it reads N c 



and 



-u 



„W0 

-L 



N c . Using this, we can estimate that the 



radius of our dominant ideal toroid behaves 
L 4tt/ 



2nN r 



WL 



(47) 



This result is, however, not directly applicable to the 
physical systems, because our model has the zero thick- 
ness of the chain. That is, every chain segment interacts 
equally with all the other segments accumulated on the 
same arc of the toroid. 

Therefore in this section, we first introduce a finite 
size effect into our Hamiltonian. We then estimate the 
mean radius of toroid and compare resulting analytical 
expression with the experiments of DNA condensation 
0,0]. Also, the mean radius of the toroid cross section 
is calculated in the end. 



A. Finite size effect 

The finite size effect of the toroid cross section can be 
approximated by the hexagonally arranged DNA chains 
with van der Waals type interactions, i.e., with the effec- 
tive nearest neighbour interactions. Namely, if the chains 
are packed in a complete hexagonal cross section, the 
winding numbers are N = 7, 19, 37, and so on. In such 
cases, the number of van der Waals interactions between 
segments can be counted by the links between neighbour- 
ing pairs in the hexagonal cross section. We then obtain 
the number as a discrete function Vdiscrete(N) . We can 
approximate it or analytically continue to the following 
analytic function V(N): 



V(N) = 3N - 2VsJ N - ~. 



(48) 



Thereupon, the attractive energy can be expressed by 
those of N loops and the rest of the chain: 



ds Vat(s)-- 



2ttW 



a 



V{N)-W\L 



2nN 



Gap(N)(49) 



where Gap(N) = V(N + 1) - V(N) is additionally in- 
troduced in order to compensate the continuity of the 
potential as a function of a. Note that, up to N = 3, we 
need not to introduce this finite size effect, since there 
is no difference between the ideal toroid and the hexag- 
onally arranged case: the number of links are the same 
in both cases. Therefore, we assume N(L) > 4 for this 
effect. Note also that the entanglement (knotting) effect 
of the chain arrangement is neglected. 



Finally, substituting ea. (|49|l into Tl[u, W], the finite 
size effect leads to the modified Hamiltonian: 



H(a) = ^ - 
2ttW 



V(N(L)) 

- N(L)) Gap(N(L)). (50) 



2ttW 
a 

aL 
2tt ' 



Most physical observables for tightly packed toroids can 
be quite accurately estimated using this Hamiltonian. 
For example, in principle, we can derive the exact value 
of the radius of the stable toroid. In fact, by the same 
analysis presented in the previous section, we obtain the 
following "asymptotic" relation of N c of the dominant 
toroid for large c: 



(2V3 cj 



(51) 



B. Mapping onto experimental data 



By r c = 



_ L 



2trN c 



, we now estimate the mean radius of 



the toroid (i.e. the average of inner and outer radii) in a 
physical system. A coupling constant of can be given 

by W = -f-^-^j^j wnere & i s the number of the electric 
dipoles in a monomer segment, each of which creates van 
der Waals interaction of the magnitude e. l m denotes the 
length of the monomer along the chain contour, taken to 
be a half of pitch per turn (of helix) l m ~ 5 bp = 1.66 nm 
in the end. Note we assume l m ~ lb- Substituting N c of 
the dominant toroid state l|51|l and the above, we obtain 



(«»)" 



knT 



(52) 



The scaling property of the first equality matches with 
the one in |25|]. Note that a coupling constant of sec- 
ond nearest neighbour is vanishingly small W2 — 2~ 6 W = 
gjW, so that one may neglect it in this van der Waals 
regime. 

We estimate the mean toroidal radius of T4 DNA in 
low ionic conditions reported in [13l |. Using L = 57/im, 
Z~50~60nm, and l m , the mean radius of the toroid is 



29.09B" 



31.29B" 



(53) 



where B = This is in good agreement with the 

experiment r c ~ 28.5 nm for 1.15. 

The same argument for the toroid formed by Sperm 
DNA packaged by protamines (L = 20.4/zm), gives 
the analytic value 



23.69B" 



25.48B" 



(54) 



which also agrees with an experimental result r c ~26.25 
nm for J5~0.85. Note that the former toroid is densely 
packed [13| , hence hexagonal assumption could be a good 
approximation. It is therefore well expected to have the 
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stronger attraction compared to thermal fluctuations, 
B > 1. The latter has a larger diameter of the effec- 
tive segment. Thus, it may well be expected to have the 
weaker but large enough interaction with smaller B to 
maintain a toroidal conformation. 

Similarly, the mean radius of the toroid cross section 
can be calculated for the complete hexagonal cross sec- 
tion with a side of n monomers: 



>/3 



1 



(55) 



where Id is the diameter of the segment. Substituting 
our relation of n with the winding number TV: n = i + 



"-/j y N — \, eq.JSSJ can be rewritten as 



2\/3 + 3 
12 



N2 



1 

8iV 



O 



N 2 



Id- (56) 



Using this general result, we now specify the type of 
interactions and derive the expression for the mean radius 
of the toroid cross section. In the case of van der Waals 
interactions, N c is given by ea. (|51|l . and therefore we have 



3V3 + 6 
12 

3^+6 
12 



(6tt)-^^ 



ke 



Id- (57) 



Note that the scaling property r cross ~ Ls is in agree- 
ment with the one in [25j obtained in the asymptotic 
limit. 

Similarly, we can formally consider the case of ideal 
toroid, although it has zero thickness. In this case with 
N c ~ c, we have 



2V3 + 3 /W 
24tt V 2 1 



2V3 + 3 T . , ,1/ ke V, 



(58) 



It should be noted that, with van der Waals interac- 
tions, the mean radius of the toroid scales as L~s while 

2 

that of the toroid cross section scales as L s . In the case 
of ideal toroid, the mean radius of the toroid scales as 
L _1 while that of the cross section scales as L. 



VII. DISCUSSION AND CONCLUDING 
REMARKS 

We have shown that the (meta-)stable toroid states ap- 
pear as the classical solutions of the low energy effective 
theory of a semiflexible homopolymer chain — the non- 
linear sigma model on a line segment. We have shown in 
this paper the complete proof of the statement that our 



classical solutions represent the general solution. The 
novelty of the model comes from the fact that the dif- 
ficulties of the local inextensibility constraint \u\ 2 = 1 
and the delta-function potential are resolved explicitly 
with our solutions in the path integral formulation (hence 
our theory goes further beyond the Gaussian approxima- 
tion). Together with our microscopic Hamiltonian, they 
lead us to the profound analytic curve of the energy lev- 
els (Fig|3J. It is also of interests that the balance be- 
tween the bending energy and the attractive potential 
creates multiple local minima at a = a c (N) for each N 
satisfying ea. (|3*2")> . One can read off from eq.{32J and 
Fig^]that the number of local minima is basically three 
or four in the case of our ideal toroid. Another work 
on multiple local minima will be mentioned shortly. In 
search of the ground state, we found that the phase tran- 
sitions occur at c = I and at c = ||, and discovered 
that such configurational transitions are governed by the 
conformation parameter c = -Sr (tt-) • The critical point 
c = H indicates the whip-toroid or whip-dominant to 
toroid-dominant transition of the first order. 

In section V-A, we have shown the stability of the 
toroid states and the validity of such phase transitions. 
We also calculated the potential barrier of the whip- 
toroid transition (rod-toroid transition), and indicated 
that, for the chains of L < 1 /im, the rod (a = 0) to 
N = 1 (meta-)stable toroid state transition is fairly un- 
likely. This would explain why it is difficult to observe 
short DNA toroids in experiments [5JJ. In section V-B, 
we have constructed the effective Green function from 
those of the whip and toroid states using perturbation 
theory at low energy. It naturally contains multi-tori 
and 'tadpole' conformations. 

Wc finally introduced the hexagonal approximation to 
count the finite size effect of the toroid cross section and 
established the mapping onto the experimental data of 
the DNA toroid radii. Our result is even quantitatively 
in good agreement with the experiments 0,|u|- Hence, 
we conclude that our theory is certainly an analytic the- 
ory of DNA condensation and of toroidal condensation 
of many other semiflexible polymer chains with effec- 
tive van der Waals attraction, or equivalently with ef- 
fective short-range dominant attraction. Here, we pre- 
sented only the comparisons with DNA condensations, 
but simply by varying the parameters l,W,L, and l m , 
our theory and results should fit to the same problems in 
similar biochemical objects. 

In analogy to the classical limit * — ► 00 in quantum 
mechanics, it is assumed that the persistence length I 
is large enough for our low energy theory to be valid. 
The local inextensibility constraint is originally given by 
some bond potential such as \(\u\ — l) 2 with the spring 
constant A. Note the low energy theory becomes invalid 
when the constraint does not hold in the Hamiltonian. 
That is, I and A should be sufficiently larger than O(l) 
so that our theory remains valid. On the other hand, as 
I approaches 1 or 0, one may be able to see the transition 
from whip-toroid phase to coil-globule or coil-rod phase. 
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This is beyond our scope in this paper, but the issue will 
be discussed at the end of this section. 

Of particular interest on multiple local minima is the 
work by Kuznetsov and Timoshenko j3Sj . Using Gaus- 
sian variational method and off-lattice Monte Carlo sim- 
ulation, they numerically calculated the phase diagram 
for a ring scmiflcxible chain in various solvent conditions. 
Note that their model Hamiltonian is the same as ours 
except for the harmonic spring term connecting adjacent 
beads. For a given stiffness, upon increasing the mag- 
nitude of two-body attraction, they found toroids with 
larger winding numbers become more stable. This is con- 
sistent with our findings as our conformational parame- 
ter c is a function of the magnitude of the attraction W . 
They also found that several distinct toroid states (i.e. 
multiple minima) can exist, which are characterized by 
winding numbers and are separated by first order phase 
transition lines (see Fig. 4 in 38]). Although they consid- 
ered a ring polymer, these facts are in good agreement 
with our analytic findings for an open semiflexible chain. 
It would be of great interest to numerically check the ex- 
istence of multiple minima for an open semiflexible chain. 

So far, we have only dealt with the whip and toroid 
conformations, but there are other configurations to be 
explored. Numerical simulations showed that a semiflexi- 
ble chain takes toroid, collapsed rod or racquet conforma- 
tions, depending on chain length, stiffness, magnitude of 
interactions, temperature, and other variables. Of par- 
ticular interest is the works by Noguchi et al. [361 137| 
and Stukan et al. [2^, who studied the dependence 
of stiffness on conformational properties. They observed 
both toroid and collapsed rod states for some intermedi- 
ate stiffness. Upon increasing stiffness they found toroid 
states are more probable. 

Collapsed rods have not been present in our model. 
One of the reasons is they are not classical solutions, 
which can only survive and become the only candidates 
for the ground state at large I or at small T. Another rea- 
son may be because the inextensibility constraint \u\ 2 = 1 
is quite strong or because collapsed rods are energetically 
less favoured. That is, our model with the constraint is 
simply in the quite stiff regime where the collapsed rods 
are less likely. Moreover, discrete nature of the chain, 
which are present in most numerical models, might al- 
low sudden hairpins although they are highly disfavoured 
in some continuum models. Indeed, when they increase 
stiffness, toroid states are more probable [371 l4lj. This 
competition in the intermediate stiffness remains an in- 
teresting open question. 

In addition to toroids and collapsed rods, tadpole like 
conformations (i.e. a toroid head with long tail) have 
been observed in the experiment by Noguchi et. al. [3?| . 
They also performed the Monte Carlo simulations and 
found that this tadpole like structure was realised only 
twice in a hundred runs. We have also dealt with a 'tad- 
pole' conformation in the effective Green function, but it 
only includes the simplest tadpole shape: a toroid with 
a single non-interacting whip. Therefore, we have not 



counted tadpoles such as a toroid with two whips at- 
tracting each other or a toroid with a collapsed rod or a 
toroid with two collapsed rods attracting each other. To 
compare them, one should first estimate the entropy of 
the whip and compare it to the toroids. Also, relative 
energy levels and statistical weights of these conforma- 
tions (toroid, collapsed rod, tadpole) includi ng repo rted 
racquet states of metastable intermediate j2^,|23,^2 will 
be studied in the forthcoming work. Note that the whip 
is defined by the elongated state at low energy whose up- 
per bound is given tentatively by E < ^j-*-. Precisely 
speaking, it should be the lowest bound for a chain to 
form a loop, which is to be explored in detail as well as 
the above. 

As for the scaling property r c ~ L v of the toroid ra- 
dius, the exponents v predicted in the literature in the 
asymptotic limit are v = | in most cases [2^, [2^, [23, OSI ■ 
This agrees with our precise asymptotic result l|52[) where 
both the parameter c and the winding number N c of a 
dominant toroid are large enough. Note our model has 
robustness in that it can treat chains of any finite length: 
a real chain is a finite system. 

However, the exponents are inconsistent with the ex- 
perimentally well known observation that the radius is 
independent of the chain length 0, 0, ^| . This might 
suggest that the real interaction is not necessarily van der 
Waals like, or at least is not a single van der Waals type 
interaction. It should be noted here that combinations of 
our ideal toroid and its finite size effect can give a range 
of v = — 1 ~ i in some region. 

Another interesting remark is that when we apply 
Coulomb like interactions to our approximation, we ob- 
serve the radius remains nearly constant as L changes. 
This will be presented in the near future [5^ . 

Finally, we remark that our model can be regarded as 
the linear sigma model at low energy where it actually 
reduces to the nonlinear sigma model. The linear sigma 
model is one of the most suitable models to describe 
phase transitions in quantum field theory. Although we 
are not formulating quantum field theory, the model ac- 
tually involves a phase transition from a constrained to 
non-constrained system, that is, from constant-length 
bonds to spring-like bonds (Gaussian flexible chain). 
This is an interesting model to be studied. Although 
we have listed some questions to solve, there are obvi- 
ously a lot of problems to be investigated. Our theory 
could also be extended and applied to the challenging 
interdisciplinary problems such as protein folding. 



Acknowledgments 

We are grateful to K. Binder, V. A. Ivanov, W. Paul, 
and N. Uchida for their stimulating discussions. Y.I. 
is grateful to K. Nagayama for his discussions and en- 
couragement, and to T. Araki for his comment. N.K. is 



13 



grateful to H. Noguchi, A. Cavallo and T. Kawakatsu for 
stimulating discussions, and to T. A. Vilgis for his ear- 
lier discussions of the field theory of globule-toroid tran- 
sition which led us to this direction. Y.I. acknowledges 
the Yukawa Institute for Theoretical Physics where this 
work is partially done during the YITP-W-05-04 Work- 
shop "Soft Matter as Structured Materials". N.K. ac- 
knowledges the Deutsche Forschungsgemeinschaft for fi- 
nancial support. 



APPENDIX A: THE DELTA FUNCTION 
POTENTIAL 

Our delta function potential expressed in the body is 

V AT {s) = -W [ da' 5 {r{s) - r{s')) . (Al) 
Jo 

This is, however, rather schematic. The precise defini- 
tion will be given below, sorting out two ambiguities in 
ea. (|Al(l . One is the definition of the delta function and 
the other is its integration contour concerning the self- 
interaction contribution. The exact form of the function 
is given by renormalising the coupling constant, the mea- 
sure, or the delta function itself appropriately, and by 
expecting some ultraviolet (short-range) cutoff e > in 
the integration contour: 



Vat (3) 



-W [ds'} 
Jo 



5{\r(s) -f(s')|) fore<s,(A2) 



otherwise it vanishes. 

First, we should interpret that the delta function is 
not three-dimensional but one-dimensional, changing its 
argument from (r(s) — r(s')) to \f(s) — f(s')|. When the 
argument of the delta function has some zeros, it has the 
following property. Say that the argument is given by a 
function g(x) then 



S(g(x)) = J2 



8(x - Xj) 
\9'{xi)\ 



(A3) 



where {xi} is the roots of g(x). Accordingly, at around 
every zero, the integration over x gives a \g'(xi) | _1 contri- 
bution. By definition, our potential should not have such 
a contribution, so that the amplitude must be normalised 
appropriately. For example, such a renormalisation of the 
measure can be achieved by 



[dx] re S (g(x)) = 



5(x - Xi) 



[°° ^ 5{x - 

/ dx S(x — Xi) 
Number of zeros. (A4) 



As given in the second line of the above expression, the 
measure is renormalised such that it cancels all the de- 
nominators of the delta function potential. 

Another example is to renormalise the delta function 
as follows: 



ds' S re (g(s')) 



dg(s')S(g(s')), (A5) 



where C is given by the one-parameter contour of g(s') 
from s' = to s' = s. Therefore, we implicitly include 
one of these renormalisations in ea. (|Al|) . 

The second ambiguity is that, when the integration 
contour is from to s, it arises a self-interaction between 
a point at s and an adjacent point at (s — e) with some 
infinitesimal positive parameter e. In order to avoid such 
a self-interaction, we have implicitly introduced an ultra- 
violet cutoff e (e > 0) in the form of Vat(s)- 



ds'S{---), 



Note that Vat(s < e) is defined as nil. 



(A6) 



APPENDIX B: SO(3) TRANSFORMATION 

The dimension of the generators of SO(3) are three: 
T l for i = 1 to 3 and global SO(3) transformation can 
be given by its exponential mapping: e 9iT where are 
some arbitrary parameters. The matrix form of the gen- 
erators on the fundamental representations are antisym- 
metric {(T l )j k } where j and k are matrix suffices run- 
ning from 1 to 3. Its adjoint representation is given by 
{T l )jk — tijk where ey& is the complete antisymmetric 
tensor. 

The SO(3) infinitesimal transformation of the bond 
vector u is given by 



Suj = g{T] k u k , 



(Bl) 



where i, k are summed over and gi are infinitesimal pa- 
rameters in this case. Clearly, the action is invariant 
under such transformations: 

(u i y(u i y~u i u l = (u l + Su^m + Sui) - \u\ 2 

= g k ((u t T k u 3 ) T + ( Ul T k u 3 )) 

= 0, (B2) 

since T k is antisymmetric (T k ) T = —(T k ). Note that 
Ui and dui have the same property under SO(3). For 
simplicity, let us take its adjoint representation Tj fe = 
Cijk and write down the transformation law in the polar 
coordinates 



)Uj = gtSijkUk, 



(B3) 
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where i, k are summed over. Substituting the polar de- 
composition with the constraint \u\ 2 = 1, one obtains 

5u x = S(sm8 u cos(f) u ) = {g 3 sm9 u sm<p u - g 2 cos9 u ) 
=> 59 u cos 9 U cos 4> u — 6(j> u sin 9 U sin <fi u 
= {g 3 sin 6 U sin (f> u - g 2 cos 9 U ) , 
5u y = (5(sin 9 U sin 4> u ) = (ffi cos 0„ — 53 sin 9 U cos </>„) 
=>■ 89 u cos sin + 5<j> u sin 0„ cos <j> u 
= {gi cos 0„ - g 3 sin U cos <j> u ) , 
Su z = 8(cos6 u ) =sm6 u (g 2 cos(j> u - gism(f) u ) 

=> S6 U = (51 sin 4> u - g 2 cos 4> u ) . (B4) 

Thus, the transformations of 9 U is shown in the last line. 



From the first and second transformations, one finds 

{second) cos <fi u — {first) sin <p u 
=> 5(j> u sin 9 U = {gi cos 9 U - g 3 sin 8 U cos <p u ) cos <p u 
- {g 3 sin 9 U sin U - g 2 cos U ) sin <p u 
8<j) u = (cot U (51 cos 4> u + g 2 sin <j> u ) - g 3 ) . (B5) 

Therefore, the SO{3) transformations are given by 

S9 U = g 1 sm(f> u - g 2 cos(f> U} 

5<t> u = cot6 u {gicos<f) u + g 2 sm<j> u ) - g 3 , (B6) 

where gi are arbitrary infinitesimal parameters which 
represent rotations around z-axis, i.e., x-, y-, and z-axes. 
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